Vector field models of modified gravity and the dark sector 
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We present a comprehensive investigation of cosmological constraints on the class of vector field 
formulations of modified gravity called Generalized Einstein- Aether models. Using linear perturba- 
tion theory we generate cosmic microwave background and large-scale structure spectra for general 
parameters of the theory, and then constrain them in various ways. We investigate two parameter 
regimes: a dark-matter candidate where the vector field sources structure formation, and a dark- 
energy candidate where it causes late-time acceleration. We find that the dark matter candidate 
does not fit the data, and identify five physical problems that can restrict this and other theories 
of dark matter. The dark energy candidate does fit the data, and we constrain its fundamental pa- 
rameters; most notably we find that the theory's kinetic index parameter n ac can differ significantly 
from its ACDM value. 



I. INTRODUCTION 

Over the past few years, a new suite of models for the 
dark sector has been proposed. They invoke a vector 
field which is normally constrained to lie along the time- 
like direction and may lead to modifications to the grav- 
itational sector [lJ-Q - Sometimes called Einstein- Aether 
models, they tend to entangle two of the main paradigms 
currently being considered: on the one hand modified 
gravity and on the other dark matter and energy. 

Vector field models are attractive because they seem 
to be able to resolve the problem of the dark sector (i.e. 
dark matter and energy) in a unified way. Most of the 
emphasis has gone into constraining vector field models 
that lead to accelerated expansion [2|-|4j although there is 
a fair amount of work for which the the vector field leads 
to a relativistic version of Modified Newtonian Dynam- 
ics (MOND) @ or can even play the role of dark matter 
[3, @. In fact vector field models seem to incorporate 
what seems to be a generic feature of relativistic modi- 
fied gravity models [7| : that it is impossible to construct 
relativistic models that just modify the gravitational sec- 
tor without introducing new degrees of freedom, which 
can then behave like either dark matter or dark e nerg y 
(although for other approaches see, for example, 

There has been significant progress in trying to con- 
strain these models. For example, at a fundamental level 
it has been shown that a broad class will lead to insta- 
bilities and the formation of caustics, signaling a break 
down of the fundamental theory [l4j]- It has also been 
shown that for a general choice of kinetic terms, these 
theories will be plagued by ghosts or tachyons p^4l7| . 
These pathologies are worrying but do not entirely rule 
out vector field models- it has been shown that modifi- 
cations to the kinetic term, for example, can cure them. 

Substantial work has been done on understanding how 
these fields in these theories behave on macroscopic 
scales, either through their interaction with matter to 
form galaxies and clusters [Isj . or on the largest scales, 



affecting the growth of structure and its effect on the 
CMB Indeed for a particular, "vanilla" version 

of the vector-field model, detailed and definitive con- 
straints have been placed on the various coupling con- 
stants!^, H3 . 

In [Hi, it was found that one of the key effects that 
vectors would have would be to modify the growth rate 
of structure. This is not surprising- theories that modify 
gravity tend to have this effect. We also found that it lead 
to a mismatch between the t wo g ravitational potentials a 



potentially observable effect [22|. In this paper we wish 
to pursue this analysis and quantify how strong these 
effects are. Although we focus on a particular (albeit 
broad) class of theories, we are interested in extracting 
general lessons from these models. We believe that much 
of what we learn by looking at these models will shed 
light on other models of modified gravity (such as, for 
example, f(R) theories [23[ and bimetric theories (24|). 

The structure of this paper is as follows. In Section 
HT1 we lay out the essential ingredients for a reasonably 
broad class of vector-like models and its background evo- 
lution, and specialize to the form used in the remainder 
of the paper. In section IIII Al we lay out the equations 
of the perturbed theory, and how we implement them 
with theoretical constraints. In section ITVl wc discuss the 
problems with modelling dark matter with the theory. In 
section [V] we find and constrain the parameters which let 
the theory act as dark energy. In section |VT1 we conclude 
and draw more general lessons about the dark sector. 



II. THE THEORY 

A. Theory Definition 

A general action for a vector field A a coupled to gravity 
can be written in the form: 



S = J d i Xy/—g 
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+ C(g ab ,A b ) 
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where g ab is the metric, R the Ricci scalar of that metric, 
Sm the matter action and L is constructed to be generally 
covariant and local. By construction Sm only couples to 
the metric, g a b and not to A a . 

We will restrict ourselves to consider a Lagrangian that 
only depends on covariant derivatives of A and we will 
consider a unit time- like A a . Such a theory can be writ- 
ten in the form: 



These equations allow us to study a general theory of 
the form presented in equation Q] with a unit time-like 
vector field. For our particular, restricted choice of K 
we have Y ab = - Cl [{\7 c A a (\7 c A b ) - {\7 a A c )(\7 b A c )} and 
Vb = 0. 



B. Background Cosmology 
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where a are dimensionless constants and M has the di- 
mension of mass. We have removed an additional C4 'ac- 
celeration' term in accordance with the transformation 
described in J25| . 

As was the case with TeVeS, the constant G may 
be a different number from the locally measured value 
of Newton's gravitational constant. A is a non- 
dynamical Lagrange-multiplicr field with dimensions of 
mass-squared. 

The gravitational field equations for this theory, ob- 
tained by varying the action with respect to g ab (see [H) 
are 



Gab — Fa. 
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where the stress-energy tensor for the vector field is given 
by 



T a b 

F K 
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Brackets around indices denote svmmetrization [48j and 
Y ab is the functional derivative 



J(K cd ef ) 
V c A e V d A f — — -4 1 
Sg ab 
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The equations of motion for the vector field, obtained 
by varying with respect to A b are 

Va(F K J a b ) + F K yb = 2XA b (10) 
where once again we define the functional derivative 



Vb 



= V c A e V d A f 



5A b 



(11) 



Finally, variations of the action with respect to A will 
fix A b A b = — 1. By inspection, contracting both sides of 
(|TU|) with A b leads to a solution for A in terms of the the 
vector field and its covariant derivatives. 



In this paper we will restrict ourselves to background 
cosmologies where the spacetime is of the spatially flat 
Friedmann-Robertson- Walker (henceforth FRW) form: 



g ab dx a dx b = -di 2 + a(t) 2 S i jdx i dx j . 



(12) 



The energy momentum content of each matter fluid (i) 
in the background is taken to be of the usual form 



T i ab = (P i +p i )uiui + P i g ab . 



(13) 



Where p 1 and P 1 are the co-moving density and pres- 
sure of the ith fluid respectively. We assume that all 
fluids have identical four velocities u a = (1,0,0,0). Fur- 
thermore we define p = ■ 

In spacetimes with FRW symmetries, the vector field 
must align with the direction dt and so the vector field 
is entirely fixed to have components (1, 0, 0, 0) in the co- 
ordinate system p^|) . Explicitly, the background value 
of the scalar K is given by: 



K 



FRW 



aH 2 
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Where H = d t ]n(a(t)) and: 

a = ci + 3c 2 + c 3 . (15) 
The Friedmann equation then takes the form fl9| : 



F 



dK \ K 1 / 2 



H 2 = 
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where p still includes only the matter components. 

We note that the combination of (Tro]) and the aether 
stress energy tensor allows us to write down an effective 
energy density p ao and pressure P ae of the aether. From 
this we may define the fractional energy density f2 ae = 
SirGpuc/SHj 2 , and the aether equation of state parameter 

^ac = -^ae / Pae ■ 
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FIG. 1: A schematic representation of the types of the late- 
time background evolution permitted by the functional form 
F = y(—K) n as a function of (71,7) for n < 1. 



The effect of the vector field on the background expan- 
sion may be see as an expansion rate dependent modi- 
fication to Newton's constant i.e. schematically (|T6|) is 
an equation of the form 3H 2 = 8nG(H 2 )p. It was found 
in that various forms of the function F permitted a 
wide variety of cosmological dynamics: the presence of 
the vector field variously leading to accelerated expan- 
sion, slowed expansion, rescaling of G, and recollapse as 
summarized in Figure [I] 



C. The functional form of F(K) in cosmology 

We must now specify the form of the function F(K) in 
equation (|6]). There is an obvious set of candidates here 
- we could attempt to be directly consistent with MOND 
and use the same branch of F as it uses on small scales. 
In appendix [X] we show that doing so would make it 
impossible to consistently generate late-time acceleration 
behavior in the background cosmology. Instead, we will 
use a simple and reasonable ansatz that works for the 
regime \K\ ^> 1 that we consider here. 

Existing functional forms for F(K) in the MOND 
regime typically are dominated for a single monomial 
term for \K\ ^> 1 (sec for instance [26j]) and so it seems 
reasonable to restrict the function to take this form: 



F = 7(-A')"" K < 
F = 7(+A')" ao K > 



where n ao < 1. This form has sufficient power to express 
a wide variety of behavior, and the parameters 7 and n 
shall be central to our further analysis. 



III. LINEAR PERTURBATION THEORY 
A. Formalism &z Theory 

We have seen that the vector field can have a signif- 
icant effect on the quasistatic, weak field limit and the 
background cosmological geometry. Significant evidence 
for non-baryonic mass persists on the largest cosmologi- 
cal scales [27| therefore it is vital that a relativistic the- 
ory of MOND can account for this. As mentioned, it 
has been argued [28[ that even in the quasistatic, weak 
field limit, a spatial tilt to the vector field may produce 
significant deviations under some circumstances from the 
local MOND force law. Similarly it was shown [ll| that 
in the context of linear cosmological perturbations the 
energy density associated with the projection of the vec- 
tor field onto surfaces of constant conformal time could, 
to a degree, act as a cosmological 'dark matter' . In this 
paper we shall comprehensively address the question of 
whether the Lagrangian © represents a viable model of 
the dark sector in light of precision cosmology. 

In the main body of the paper we will consider scalar 
perturbations. In Appendix [B] we derive the equations 
of motion for the vector field's two divergenceless vector 
modes as well as the gravitational wave tensor modes (the 
speed of propagation of which is in general modified by 
the vector field). Requiring the stability of these modes 
puts constraints on © (see Section fill Bp . 

We shall work in the synchronous gauge (see for in- 
stance [29| ) and so the metric takes the following form: 



g flu dx IJ, dx v = -a 2 dr 2 + a 2 [7^ + h,, d.r'd.r' (20) 

where r is conformal time, 7^ is a spatially flat space- 
like 3-metric perturbed by hij which is built from two 
scalar potentials 77 and h: 



M*.t) = / d 3 ke^ x [kk J h(k,T)+(k l k 3 --5 lJ )6i 1 (k,T)} 

(21) 

Similarly we will expand the aether field as: 



A" = -(l,diV) 
a 



(22) 



(19) 



The zeroth component of the aether field is, by virtue 
of the gauge choice and the constraint, fixed as equal to 
1 up to second order in perturbations. 

A slight complication in the field equations arises be- 
cause of the presence of the function F which depends 
nonlincarly on the scalar K. We assume that for modes 
of interest one may consistently regard the perturbation 
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to K as being much less than unity. Thus one can expand 
K as K = K° + K e and F as F = F° + F° K K e where 
K e <C 1 where the superscript denotes the quantity cor- 
responding to a function's background value. Explicitly 
we have that: 



K° 



, an 2 

'a 2 AP 

2 h' 
-±K\k 2 V--) 



(23) 
(24) 



Where primes denote derivatives with respect to con- 
formal time and we have used the conformal Hubble pa- 
rameter H = H/a. We also make use of the following 
identity: 



(K°)"H = -4K Q (H 2 -- — ) 
2 a 



(25) 



Henceforth we will drop the superscripts on K° and F° 
i.e. K and F shall be assumed to represent background 
values of the fields. 

Towards simplifying the form of the equations, we will 
will rather use the field £ instead of V, where £ is defined 
as: 



(26) 



For further compactness of notation we define the vari- 
ables: 



l + 2^if| a 



K 



1 + 



F; 



K 



(27) 
(28) 



The vector field equation of motion (fTU| becomes 



= c l {l + c l3 F K )^P^ + 2Hc l {l + c l3 F K )^ F 



F K 



Fk 



+[2ci(1 + c 13 F K )(— - H 2 ) + 2(ci + a)(H 2 - ™) 
a 2 a 

+c 1 c 13 {F KK K r y + + 2ci 3 )fc 2 ]£ 

+ (ci + <i)T/ + (c! + a)^(n 2 - \—){h! + 6r/) 



3d (F K z f y 

2 k 2 F K 
The relevant Einstein equations are: 
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k 4 {a + 2c 13 )F K £ (30) 



and 



(1 + -c x F K )(Hti - 2fc 2 ?/) = -8TrGa 2 5T° 



--F K ( Cl + &)6Hr)' 
-2aF KK KUk 2 t 

-^-(« 2 (i + c 13 ^)e)' 



(31) 



Where we have used the fact that: 

(h + 6r])" + 2n(h + 6r]Y -2k 2 rj = -3E / 

+2c 13 fc 2 [F K (2?^ + £') 
+F A - K A'£] (32) 

The functions 8T°j and <5T° are the first order pertur- 
bations to the corresponding components of the matter 
fields' stress energy tensors. Summation over field species 
is assumed. The field £/ is the scalar component of the 
total fluid shear i.e. £/ = —8nGa 2 (kik 3 — ^S 3 ^ 1 - and 



ST\ 



A gauge invariant formulation of 



the theory's equations may be found in [3C 



B. Parameter Constraints 

We can immediately see a number of constraints on 
the Ci and the form of the function F. From (p?5]) it can 
be shown that for in the limit of timescales shorter than 
a Hubble time the quantity: 



O (9L 

r<z — v 2 

L>S — - 



C13) 



3 Cl (l + c 13 F K ) 



(33) 



can be interpreted as the squared sound speed of the field 
£. The avoidance of exponentially growing subhorizon 
modes dictates that C| should be positive definite. 

Similarly, one may consider the field equations of the 
two divergenceless 'vector' modes of the vector field and 
the two transverse traceless 'tensor' modes of of the met- 
ric. Each respectively has a squared sound speed func- 
tion (named C v and G\ respectively) which, as in the 
scalar case, should be positive definite. These functions 
are calculated in Appendix iBl and are as follows: 



CJrp 



1 



1 + c 13 F K 



G 



F K 2 Cl +F K (c 2 -c 2 ) 
2ci 1 + c 13 F K 



(34) 



(35) 



Note that the gravitational wave tensor modes now 
^enerically have a time dependent speed of propagation. 
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Collectively the three positivity constraints imply the 
following constraints on the q and function F parameter 
space: 



(a 



l + F K c 13 > 
-2ci 3 )/2ci > 



F K (1 + F K 



r 3 
< 1 



2ci 



> 



(36) 
(37) 

(38) 



Now we turn to the Einstein equations (|3U| and (f3"TT) . 
Terms in the perturbed vector field stress energy ten- 
sor ST a b may contain terms proportional (up to a time- 
dependent function of the background fields) to SG a b and 
indeed this can be seen in (|30|) and (|3T|) through the ap- 
pearance of terms proportional to the Cj on the left hand 
sides of the equations. If the effect £ is negligible then 
the effect of the aether stress energy terms proportional 
to components of the Einstein tensor may always be ab- 
sorbed into a redefinition of Newton's constant as a time- 
dependent effective gravitational coupling. Thus, even if 
£ has comparatively little effect, there may be a consid- 
erable modification to the link between the matter fields 
and the gravitational field. The resulting gravitational 
couplings should be greater than or equal to zero other- 
wise the gravitational field will interpret normal matter 
as violating energy conditions, and so risking the appear- 
ance of instabilities. This restriction implies the following 
constraints: 



(1 - -aF K ) > 



(1 + ^Fk) > 



(39) 



(40) 



Throughout our analysis, we will only consider regions 
of the model's parameter space which satisfy these con- 
straints. 



C. Computation 

To study the effects of the vector field in detail we 
have modified the structure formation Boltzmann code 
CMBEASY [HI. We add a Newton-Raphson solver for 
the Hubble parameter, with added aether components. 
The perturbation evolution is also modified to include 
the aether components £ and and their contribution 
to the density, pressure and shear perturbations. We 
also include the altered metric perturbations in the cal- 
culation of the CMB source function. We use adiabatic 
initial conditions [HJ. Since Boltzmann codes are very 
highly optimized for ACDM models, care must be taken 
to ensure that modifications are performed in a consistent 
manner - for example an unmodified Fricdmann equation 
is often assumed for computational efficiency. 



To explore parameter spaces, CMBEASY is coupled 
to a Monte-Carlo Markov Chain (MCMC) engine [33| . 
We extended this engine to include our new Aether pa- 
rameters: r2 ac , c+ = (ci + c 3 ), c_ = (ci - c 3 ), C§, n ao 
and M. We include the full ranges of these parameters 
by allowing the kinetic term to take two branches, for 
positive and negative K as in (|19p. 

We constrain models using both CMB data from the 



WMAP experiment 
the SDSS survey (35 1 



34J and large-scale structure from 
though not always at the same 
time. As we shall see, some Aether models are ex- 
tremely poor fits to the combined data sets; to illustrate 
such problems we want to find models that fit only the 
large scale structure data. In regimes where fits are ex- 
tremely poor, MCMC does not work particularly well. 
We ameliorate such situations by running larger num- 
bers of shorter Markov chains and re-starting from their 
best-fit positions, and sometimes by abandoning MCMC 
altogether and simply performing random searches for 
good parameter combinations. 



IV. THE VECTOR FIELD AS DARK MATTER 
A. Background evolution &z doppler peak positions 

An unusual property of the model considered here is 
that at the level of cosmological perturbations the field 
can mimic a perturbed pressureless fluid in the formation 
of large scale structure whilst behaving entirely differ- 
ently in the cosmological background jl9j . In this paper 
we would like to consider not just large scale structure 
but also other cosmological probes. The anisotropies in 
the CMB temperature are sensitive to the background 
dynamics and perturbed dynamics of source of cosmic 
mass discrepancies in a largely distinguishable manner. 

If the aether plays the role of dark matter only in the 
perturbations, then the background expansion is dark 
energy-dominated at an earlier time. For a given Hq 
this reduces the expansion rate of the universe between 
recombination and now and so decreases the angle sub- 
tended on the sky by given distance at last scattering. 
This moves the CMB doppler peaks to higher I. Al- 
though the aether can give a suitable time dependence 
to the effective gravitational coupling G in the Fried- 
mann equations so as to yield the same expansion rate 
as ACDM, the necessary functional forms of F are ex- 
tremely contrived. For instance these forms essentially 
contain a new constant scale roughly equal to the Hub- 
ble parameter at matter-radiation equality. Even in the 
event of such a construction, it may not be possible for 
the squared sound speed of vector field perturbations Cg 
to remain sufficiently small as to behave like cold dark 
matter in perturbations 

The other input to the peak position is the sound hori- 
zon at last scattering. The physics of this is sufficiently 
robust that changes compensatory to the alteration in the 
distance to last scattering arc not feasible without excep- 
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tional fine tuning of the perturbational behavior of the 
aether. Indeed, it has been argued that only additional 
components that behave like non-relativistic matter in 
the background might fix this problem (36j . 

Could there indeed be such a non-relativistic matter 
in the background which allows acoustic peak positions 
which are consistent with the data whilst leaving the 
aether in its role as the seed of the formation of large 
scale structure? A appropriate candidate would appear 
to be massive neutrinos. A suitable mass of such particles 
so as to account for the right effective contribution to the 
dust component of the background typically implies that 
the neutrinos themselves are unable to clump on small 
enough scales so as to be a good candidate for all of the 
dark matter. Such a solution as has been proposed in 
[33 and M. 



B. Perturbation evolution 

We now argue that even with the inclusion of massive 
neutrinos, one generically expects the aether to have an 
unacceptable influence on the large scale CMB anisotropy 
if it is to also play the dominant role in structure forma- 
tion. 

The first requirement for successful perturbation evo- 
lution is that structure can form at all. One necessary 
condition for this is that the sound speed of the structure 
seed not be too large, since this would wash out struc- 
ture. We require that the sound horizon in the model 
be less than the smallest scales where structure can form 
linearly: Csfc max r < 1, where fc max ~ 0.2/i/Mpc. For 
matter power observations at r ~ 3 x 10 4 , the present 
epoch, this yields C5 < 10~ 4 . 

There are two underlying physical processes that fur- 
ther constrain the models. 

The first is a change to the growth rate of perturbation 
amplitude. This can cause discrepancies between the am- 
plitudes we expect in the matter power spectrum and the 
CMB, since the evolution between the two is different. It 
can also lead to an integrated Sachs- Wolfe (ISW - see 
below) effect during the matter era since <J> will accrue a 
time dependence. 

The second is the increased presence of a $ — "J metric 
shear. This also leads (directly) to a matter era ISW. 



where r* is the conformal time of last scattering, tq is the 
conformal time today, d« is the comoving radius of the 
surface of last scattering, and $ and $ are the conformal 
Newtonian gauge gravitational potentials. The integral 
in flU) is the integrated Sachs- Wolfe (ISW) effect. Writ- 
ing the integrand of (gTJ as + = — + 2$', 
we can see this as time derivatives of a shear part and a 
growth rate part. 

In the standard cosmological model, the field $ has 
negligible time dependence during matter domination. It 
gains a time dependence only when the background starts 
accelerating, and only then can the resulting growth rate 
ISW contribution be considerable. In the model consid- 
ered here, the situation may be rather different. It was 
found in [19| that substantial contributions to the ISW 
may occur even during the matter era. For the aether 
field to seed structure formation the field £ must have a 
suitable growing mode solution in the matter era. Typ- 
ically the corresponding spatial curvature perturbation 
fc 2 <3? will then have a time dependence via the Poisson 
equation. 

The shear part can also gain a time dependence in the 
aether model, which in the ACDM is very small even 
during acceleration. 

Each of these effects depend on the functional form of 
F, the time-dependence of the £ growing mode and the 
choice of the parameters c^. It is extremely challenging 
to find combinations of the parameters which allow for 
a realistic growth of structure whilst making the ISW 
acceptably low. This is most easily illustrated by con- 
sidering the theory TeVeS [24| which has many of the 
same properties as the model considered here. It may 
be shown that TeVeS can be written as a single met- 
ric theory with a timelike vector field of unfixed norm 
J3j|,|4(j. As in the model considered here, the longitudi- 
nal component of the vector field can source the growth 
of structure j37j], dJ. We will call this field V T . We con- 
sider a matter dominated era where Vr is responsible for 
the dominant source in the Poisson and shear equations. 
These equations respectively then arej4lj|: 



fc 2 $ « -f s Hk 2 V T - ^-k 2 V+ 
fc 2 (*-$) » f s k 2 {2UV T -Vi) 



fs(r) 



A 4 



(42) 
(43) 
(44) 



Observable 1: ISW 

Under the assumption of adiabaticity we have that the 
anisotropy in the CMB, AT(n)/T on large scales in a 
given spatial direction n is given by 



where Kb is a positive constant of the action and A 4 
is the norm of the vector squared again (equal to unity in 
the fixed norm case, but in TeVeS the deviation of \A 2 \ 
from unity is essentially the background variation of the 
'scalar field' degree of freedom). Therefore: 



^ * -!*(T, )t W) - [° + #)[f, (to - f)n] * + * = "(/.W + *bM ^ 

(41) In this era the vector field equation is: 
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Vr+&i— + &2^ = (46) 

where 

61 = 2(3 - A 4 ) (47) 

6 2 = 2(2-i 4 ) + —- (l-A 4 ) (48) 

and is a source term. 

For this situation to arise, there must be a growing 
mode in Vt [III- Therefore we require that b 2 < 0. The 
function A 4 will be rather close to unity [H| so 

b 2 ~ 2 + 12 f,(T)/ K B (49) 

Therefore we require that f s {r)/KB i5 —1/6 (recall- 
ing that A's is a positive number). By ([450 we see 
that the contribution to the ISW will be proportional 
to / s (r) + Kb- Meanwhile we see from (|42p that the 
comparative contribution to Poisson's equation is also 
independently weighted by f s and Kb- If there is no 
time dependence in the Poisson equation (|4"2"| due to the 
vector field then the vector field will make no contribu- 
tion to the ISW effect. If there is a time dependence in 
the Poisson equation, there may yet be no contribution 
to the ISW if / s (t) ~ —Kb between last scattering and 
now (though the overdensities of the baryon field will 
then not generally grow as a, thus contributing to the 
ISW effect). This would be consistent with the condition 
for a growing mode in Vt but it does not guarantee that 
the resulting growing mode would be suitable. 

Indeed, it was found [37| that in seeding the growth 
of large scale structure in TeVeS there was a signifi- 
cant danger of incurring unacceptably high temperature 
anisotropics in the CMB on large scales. Although in- 
volving a larger number of terms, the same reasoning 
carries over to the model considered here - i.e. param- 
eters which allow realistic structure formation will typi- 
cally lead to an unacceptable ISW effect. This is vividly 
illustrated in Figure 0] which shows best fit models as 
compared to SDSS large scale structure data. In every 
case, the corresponding temperature anisotropy displays 
a dramatically poor fit to the data at low I. 

Observable 2: Amplitudes 

The ratio of the observed amplitudes of the CMB 
anisotropy and the matter power spectrum is consistent 
with a growth rate proportional to a (though see [43j ) . 
Any uncompensated change to this growth rate in the 
aether model over this period would lead to a different ra- 
tio. The bias parameter between the galaxy distribution 
and the underlying density field can be used to rectify this 



difference, but only if the change is relatively small and 
unphysically large bias parameters (larger than ~ 10) are 
not required. 



C. Summary 

We have seen that although the vector field may play a 
number of the roles that dark matter plays, it seemingly 
cannot do all at once [49j | . 

The position of the acoustic peaks in the CMB temper- 
ature anisotropy should be taken as a strong indication 
of additional nonbaryonic nonrelativistic mattter present 
in the universe during matter domination. Such an effect 
can be achieved in this model by a rcscaling of the value 
of Newton's constant. However, this rescaling cannot 
persist into the radiation era j44[. Thus, the functional 
form would have to be approaching the rcscaling solution 
only after one would expect matter (including cold dark 
matter)-radiation equality to happen. 

This implies the presence of a new scale in the theory, 
roughly corresponding to the Hubble parameter H eq at 
this time. It seems fair to say in general that a model 
such as that considered here, is more likely to be a cosmo- 
logically viable candidate for dark matter if the scale M 
in the theory is closer to H eq and not Hq. It is tempting 
to speculate whether a theory where the scale M itself 
is dynamical may find more success, but that will not be 
explored in this work. 

Furthermore we have seen that even if the background 
is consistent with observations, the effect on the evolu- 
tion of perturbations may be unacceptable, notably ei- 
ther through the ISW effect or comparing the respec- 
tive amplitudes of the CMB anisotropy and matter power 
spectrum today. 

D. Example Problem Spectra 

It has previously been shown that the Einstein Aether 
can produce acceptable matter power spectra with cer- 
tain parameter combinations [19[ . Here we show that 
such combinations do not provide an acceptable fit to 
CMB measurements. Despite extensive searches we have 
been unable to find any parameter set within the model 
that does fit the WMAP data well; this is entirely in line 
with the problems discussed above. 

We use a parameter set which is consistent with BBN 
limits on fish 2 and the HST key project measurement of 
Hq. The standard cosmological parameters are: Vlbh 2 = 
0.0193, n s = 0.83, H = 89.3km/s/Mpc, n c h 2 = 0. The 
new Aether variables are c+ = —4.72, c_ = —6.11, n ao = 
0.34, ft ac = 0.82, M ae = 111.3 km/s/Mpc with c 2 set by 
requiring zero sound speed. This parameter combination 
is in no sense optimal, but it does provide an illustration 
of all the problems that arise here. 

Figure 2] shows power spectra from our modified Boltz- 
mann code for this parameter set. The matter power is 
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a realistic fit to the SDSS data (this was the criterion 
for our choice of parameters). The CMB spectra shows 
various problems. In the low-£ regime a large ISW effect 
is clearly present, destroying the fit at large scales, as 
described in section (|IVB[) . The positions of the peaks 
are poorly fit by the model, as expected and discussed 
in section (jIV Aj) . Finally, in the plot we have rescaled 
the amplitude of the matter power spectrum by a factor 
0.02, corresponding to a galaxy bias of 0.14 in order to 
reconcile the relative amplitudes of the two spectra with 
the data; such a scaling is unphysically small. This cor- 
responds to the changed growth rate described in section 
PVB[) . All these effects cause severe problems when at- 
tempting to simultaneously fit the CMB and large scale 
structure. 

Figures [5] and [3] illustrate the sources of the extreme 
ISW effects shown in Figure HJ the time derivatives of 
the metric quantities plotted create an ISW effect as 
shown in equation (|4"Tj) . The onset of background ac- 
celeration in each case is marked by a turnover in the 
curves at late time. The GEA universe exhibits a dra- 
matically increased |$ — ^\ and time dependence of |$| 
during the matter era as compared to the ACDM uni- 
verse. Although the |$ — \t| has a smaller magnitude 
its time dependence can be significant for the total ISW 
effect. Note that values of r between the two universes 
do not correspond to the same physical time or redshift 
since the universes expand at different rates. 
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FIG. 2: Exotic behavior of metric potentials for k ~ 10~ 
Mpc -1 . The panels show the fields |$| and |$ — >3>| for a 
ACDM universe (green solid line) and GEA universe (dashed 
blue line) as a function of k(r — r r ) where r r is the conformal 
time of recombination. 



V. THE VECTOR FIELD AS DARK ENERGY 
A. Dark Energy Regime 

As shown in figure [1] our vector field can produce late- 
time acceleration and so play the role of dark energy. 
Indeed, for the form of the vector field used here, as the 
index n ao — > the theory becomes the same as a cos- 
mological constant for both the background and the per- 
turbations. Since the model can fit the data well we can 
use our MCMC engine to find constraints on the param- 
eters of the vector field, telling us exactly how close to 
the n ae = cosmological constant case the theory must 
be to fit the CMB and LSS data. 

If the model can fit the data only extremely close to 
«ac = then it does not provide a compelling alternative 
to the cosmological constant. If, on the other hand, there 
is significant flexibility in the model and no fine tuning, 
or if it can provide a better fit than ACDM, then it is 
somewhat more interesting. 

In this section we will consider the resulting back- 
ground evolution, CMB temperature anisotropy, and 
matter power spectrum for a universe containing the vec- 
tor field, cold dark matter, the conventional matter fields, 
and no cosmological constant. The acceleration will arise 
solely from the vector field's modification to the Fried- 
mann equation. 
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FIG. 3: Equivalent of Figure for k ~ 10 -4 Mpc' 1 
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FIG. 4: Matter power (top) and CMB (bottom) power spectra 
for the hCDM (dashed green) and typical GEA (solid red) 
models, with WMAP and SDSS constraints. 




FIG. 6: Constraints on the kinetic term power law index pa- 
rameter. 
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FIG. 5: Constraints on the three coupling terms of the theory. 



B. Constraints on aether dark energy from data 

Our MCMC generated the constraints on the vector 
field parameters shown in figures [5] to [8] these curves 
are the smoothed histograms from our combined Markov 
chains. 

The most important trend evident in these results is 
closeness of n ae to the A value of zero. We find the best 
fit value n ac = 4.2 ■ 10 -3 , with a 95% upper limit n ao < 
0.126. 

The best fit value in the MCMC run is very close to 
the ACDM likelihood of the same data, at the cost of 
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FIG. 7: Constraints on the vector field sound speed parame- 
ter. 



six extra parameters, meaning that it is unlikely to be 
favored by any model comparison exercise. It does, how- 
ever, demonstrate the validity of modified gravity-related 
dark energy candidates. 

Having obtained these constraints we can determine 
their origin. There are two ways in which the vector field 
must behave like A to provide a good fit. The first is 
that the late-time acceleration should be close to that 
given by A. The second is that any perturbations in the 
field (which are not present in A) should not affect the 
observable spectra. 



C. Constraint origins - acceleration rates 

The consistency of the acceleration of the universe with 
the cosmological constant equation of state w = — 1 is 
being measured with increasing precision in supernova 
and baryon acoustic oscillation experiments (which are 
beyond the scope of this paper). Here, they will be con- 
strained by the late-time ISW effect induced by dark en- 
ergy, and by the perturbation growth rate. 
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FIG. 8: Constraints on the parameter a = ci + 3c2 + cz 
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FIG. 9: The vector field's equation of state as a function of 
redshift z, for various values of the kinetic term index n ae . 



We can assess how closely vector-induced acceleration 
mimics A-drivcn expansion at late times with the equa- 
tion of state w ae of the vector field in the background: 

d ' 2 . f 1 

(50) 



-1 - 



1 



3# 2 A (f ) 



For the monomial form of F(K) we have that: 



= -1 - 



2n 



H 



3(2n- 1) H 2 



(51) 



Thus the equation of state will generically deviate from 
— 1 whenever n ^ and so the acceleration for these val- 
ues will not be degenerate with a cosmological constant. 
We see immediately from equation (|51[) that w ae (r) < — 1 
for < n < 1/2 and w^t) > -1 for n > 1/2. This is 
clearly visible in Figure 



D. Constraint origins - perturbation evolution 

Even if the background expansion is rather close to the 
the ACDM model, the evolution of perturbations need 



not be. This is most easily illustrated by considering the 
Poisson equation on large scales (see [l[ for a derivation) . 
On these large scales there is a time-dependent re-scaling 
of the the metric perturbation $, which we can cast as a 
modification of the effective gravitational constant G: 



b c ff = 



G 



1 



(52) 
(53) 



where we have assumed that terms proportional to the 
velocity divergence are ignorable and provisionally con- 
sidered the effect of the field £ to be subdominant. 

The Friedmann equation may be used to cast the above 
equation in a more familiar form by eliminating he back- 
ground pi in favour of background expansion rate of the 
universe and the time-dependent fractional energy den- 
sity Oi(-r) = 8ttG P (t)/(3H(t) 2 ). This yields: 



k 2 § 



,G 



(i) 



(0) 
eff 



(54) 



where 



G 



(0) 



G 



K~ 



(55) 



The n = ACDM Poisson equation may be cast in this 

G/(l - A/(3tf 2 )). 



form by taking = G and G 



,(0) 
r cfl 



For the case where n differs from 0, the function G^ will 
generically possess a time dependence during the back- 
ground evolution. Therefore the link between the time 
evolution of the functions G ®g , Si and $ will differ from 
the case where acceleration is caused by a cosmological 
constant. We may thus expect the ISW effect to be of a 
non-standard form. This is vividly illustrated in Figure 
[TU] where it can be seen that for a given set of (q, 7, M), 
variation of n results in a considerable variation in the 
large scale CMB temperature anisotropy. Also evident is 
a variation in the matter power spectrum amplitude with 
n, evident on all scales. 

Variation of parameters other than n could also have a 
significant impact on the success of the models. Given the 
results of the previous section, it seems unlikely that any 
influence of the field £ would tend to improve the mod- 
els. This indeed seems to be the case. Figure ITTI depicts 
various models where the function is varied C| for fixed 
values of the other parameters. In particular, n takes the 
value 0.1. The function C| is ultimately a measure of the 
ability of the field £ to sustain any homogeneous growing 
behavior for (fcr) > 1; higher values will tend to limit 
the effect of the vector field to larger and larger scales. 
The sets of parameters were chosen so such a growing so- 
lution indeed existed on superhorizon scales. The figure 
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FIG. 10: Matter power (top) and CMB (bottom) power spec- 
tra for the various GEA dark energy models, with WMAP 
and SDSS constraints. The power-law function's exponent n 
is varied 



FIG. 11: Matter power (top) and CMB (bottom) power spec- 
tra for the various GEA dark energy models, with WMAP 
and SDSS constraints. The squared speed of sound of vector 
field perturbations is varied. 



indicates that a growing £ field will indeed have deleteri- 
ous effects on scales where it is not suppressed. There is 
a significant ISW effect evident for the red curve CMB; 
the corresponding model must be considered as being at 
the edge of acceptability. The corresponding large scale 
matter power spectrum exhibits exotic oscillations, en- 
tirely unrelated to the baryon-acoustic oscillations which 
occur on other scales. Their presence is in this model is 
thus reflective of dynamics in the dark energy sector. 

VI. CONCLUSIONS 

A. Being Dark Matter is hard 

Generalized Einstein- Aether can, with different param- 
eter choices, resemble dark matter in some important 
ways but never all of them at once. Specifically, the new 
degrees of freedom introduced by the model may con- 
spire to identically replicate one or more but not all of 
the following properties of cold dark matter. 

1. Background dynamics 

To accomplish identical background dynamics to cold 
dark matter, one must introduce considerable fine tun- 



ing into the function F of the theory. Specifically the 
function must change form either side of (dark) matter- 
radiation equality. The parameter we tune to make this 
happen is in the action itself, unlike the usual case where 
we simply alter the abundance tt c . Changes to fit cosmo- 
logical observations can therefore have a larger impact on 
the small scale behavior of gravity. 



2. The speed of sound 

If the speed of sound is too high then structure can- 
not form on small enough scales. In this model we may 
reduce the sound speed to be close to zero, at the cost 
of one of our parameters. When designing new gravity 
theories this is perhaps the easiest structure formation 
constraint to investigate, and it should be examined to 
see if it conflicts with other constraints needed to make 
the theory useful - for an example, see [45j |. 



3. Growth rate of 'overdensity' 

Theories of modified gravity designed to replace dark 
matter must necessarily have growing modes of fluctua- 
tion in at least one of the new degrees of freedom they 
introduce, in order to sufficiently source gravitational col- 
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lapse and structure formation on scales within their own 
sound horizon. There remains some flexibility in the 
perturbation growth rate, since the bias on the galaxy 
power spectrum measurements is a free factor. As mea- 
surements of weak lensing (which samples gravitation di- 
rectly) and semi-analytic models (which predict bias) im- 
prove this freedom will be reduced. 

4- Absence of anisotropic stress and contribution the 
cosmological Poisson equation 

A sufficiently small anisotropic stress associated with 
the vector field may be implemented by fine tuning the 
parameter C13 to be very small (see equation (|32[1 ). How- 
ever, as was discussed in Subsection II V B[ even this will 
tend to come at the expense of other desired behavior of 
the field. An appreciable time variation of the anisotropic 
stress over the time from last scattering to today can re- 
sult in a very poor fit to the low-£ CMB Cg. This problem 
is likely to be common in theories of modified gravity. 

As we have seen, it is a combination of time variation 
of the anisotropic stress and time variation of the field 
<f> via the vector field's effect on the cosmological Pois- 
son equation that contribute to the ISW effect. Though 
both effects are absent in the cold dark matter case, one 
may imagine both effects being present in the vector field 
model but being of equal and opposite sign. As with 
the case of the isolated anisotropic shear contribution, it 
seems this is not possible whilst maintaining the other 
constraints like the existence of a growing mode. 

5. Effective minimal coupling to the gravitational field 

Even if the vector field growing mode gives an ap- 
propriate (dark matter-like) contribution to the Poisson 
equation, the link between the overdensity and the cor- 
responding $ can differ from the CDM case. This differ- 
ence can come from curvature terms in the vector field 
stress energy tensor, and its main consequence is a time- 
dependent rescaling of Newton's constant G. Thus $ 
may gain a time dependence during the matter era even 
there is a completely standard dark matter contribution 
to the Poisson equation. The converse may also be pos- 
sible: the time dependence of an incorrect Poisson con- 
tribution could be counteracted by a time dependence of 
the effective G. 



B. Being Dark Energy is easy 

As a model for dark energy the Generalized Einstein- 
Aether theory is more successful: we have obtained con- 
straints on its parameters and found that it generates 
spectra that fit the data across a reasonable range of its 
parameter space. It is clear that modified gravity ap- 
proaches to explaining late-time acceleration are viable 



and can provide motivated explanations for dark energy 
(though this theory retains the co-incidence problem in 
the guise of the parameter M). 

1. Closeness to A 

The most interesting constraint on this branch of the 
theory is on the parameter n ao and is shown in figure 
IH1 In some sense this parameter describes how closely 
the theory mimics A (which has n ac = 0). The fact that 
this parameter is rather free, n ao < 0.126, (95% CL), 
is consistent with the fact that a wide variety of other 
theories can also explain dark energy: present structure 
formation data is not very informative about the nature 
of dark energy, and deeper require expansion probes like 
baryon acoustic oscillation and supernovae. 

2. Sound speed 

The other notable constraint on Einstein- Aether dark 
energy, which may extend to other modified gravity ap- 
proaches, is the limits on the sound speed, illustrated in 
figure [7] As shown in figure [TTJ an incorrect sound speed 
can lead to large scale oscillations by modifying the other 
parameters of the theory and permitting a growing mode 
excitation at late time. 



3. Other constraints 

The other constraints on the theory (which are easily 
fulfilled by choosing the c parameters) come from ensur- 
ing that no growing mode can disrupt the power spectra, 
that the acceleration is close to the ACDM value, and 
that the value of effective G remains positive at all times. 

C. Future issues for modified gravity and structure 
formation 

Because ACDM is such a good fit to current cosmo- 
logical data, modified gravity will never be favoured in a 
model comparison exercise using only current data about 
linear structure. It is only in combination with physics 
on galactic and smaller scales that it can be be persua- 
sive. This work highlights a few issues for future model- 
building in this vein. 

The Generalized Einstein Aether model is a member 
of a class of models in which the scale M ~ Ho associ- 
ated with Dark Energy is visible to dark matter. Many 
of the issues raised here will be relevant to any such mod- 
els which try to use a dark matter scale consistent with 
small-scale modifications to gravity. 

A combination of probes sensitive only to the back- 
ground (like Type 1A supernovae) and to the behavior of 
cosmological perturbations is needed to fully constrain 
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these theories. For example, a value n ao = 0.3 has 
w(z) ~ 1 at low redshift but is ruled out by our con- 
straints. Similarly there are models with n ao ~ 0.75 that 
provide reasonable spectra, but they are ruled out by 
w{z) constraints. 

There are, of course, a number of extensions to the 
theory and features of it that could be change; we could, 
for example, allow M to vary dynamically, or add more 
terms to the kinetic component F in equation (|19|) . There 
are also myriad possibilities in more changing more gen- 
eral aspects of the primordial conditions or cosmologi- 
cal parameters: what happens if we add add tensors? 
Can we include an isocurvature mode? Would massive 
neutrinos help? Or curvature? This leads us to a key 
caveat that applies to this and all similar work constrain- 



ing new physics with linear structure: a simple constraint 
from the data alone is worthless, since any of the numer- 
ous other parameters we could change might conspire to 
counteract whatever problem it solves. We need a phys- 
ical explanation of a constraint's origin to understand 
whether it is robust to the cosmologist's tinkering. 
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The Friedman equation then becomes: 



H- 



ci 



H 



3 a 



1 + - - o" -ET — = &7tG P ( Ag ) 



Cl 



Appendix A: The MOND regime F(K) cannot yield 
late-time acceleration 



The MOND branch of F applies in the regime < 
if C 1; outside this range MOND does not prescribe its 
form. Inside that range the MOND value is: 



2 3 

dF = -2K+ -K* 
3 



(Al) 



(0 < K < 1) 
The weak field limit of the Einstein Aether theory is: 



K 



WF 



-Cl- 



M 2 



(A2) 



where <f> is the conformal Newtonian potential. The weak 
field MOND limit is: 



f/ MOND 
A WF 



(A3) 



Clearly we can equate these by setting M 2 = — ciOq. Us- 
ing this relation with equation (|14j) gives us this value for 
the present-day cosmological K: 



FRW^O 



(A4) 



There are now two cases: a/c\ < and a/ci > 0. In 
either case we need the measured values of H and a ; 
they are suggestively similar: 



H = Kfl 
K PS 6 



(A5) 



where 

If a/ci < then K > and we are directly in the 
MOND regime. Then we obtain the modified Friedman 
equation: 
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Cl 



8irGp (A6) 



A self-accelerating solution to this is only possible if the 
quantity in brackets is positive definite, so that H — > 
const as p — > 0. This could only happen if a/cl < — 1, 
but that would violate our requirements that if < 1. 

In the other case a/c\ > we must extrapolate the 
MOND form of F in equation (|J2]) to K < 0. To make 
the extension continuous across the K = 0, we should set 
F{-K) = —F(K), so that the MOND form becomes: 



Ci-F = -2K 
(0 < -K < 1) 



2 / r „2 



(A7) 



which docs have an accelerating solution. Unfortunately, 
solving this equation at the present day for reasonable 
values of fl m shows that it requires K ~ 2, which again 
violates our requirement that \K\ <S 1. Having exhausted 
our other options we are forced to require K > 1, outside 
the true MOND regime. This is consistent with another 
separate analysis of the solutions of (fT6| {42}. It seems 
likely that if the vector field is responsible for the late 
time acceleration then it is a result of behavior of the 
function away from the MONDian limit. There is still, 
though a role for the near numerical coincidence of a® and 
Ho - it lessens the fine tuning of the other parameters in 
F in order for the acceleration to happen at suitably late 
times. 



Appendix B: Vector And Tensor Modes 

In this appendix, we provide the perturbed equations 
of motion for the vector and tensor modes of the vari- 
ous fields appearing in these models, namely the metric, 
the vector field and the Lagrange multiplier. The latter, 
being scalar, only have a spin— mode. The other two 
fields, gij and A 1 are perturbed as 

ds 2 = —a 2 dT 2 + a 2 BidTdx l +a 2 {^ij + hij)dx l dx : ' 
hij = 2d(iEjj + 2Eij 
"1 A' 
a' a 



(-a,aVi) . 



(Bl) 



where the different fields introduced satisfy 

d i B i = = SflEij = d t E l = E\ and d i A i = 0. (B2) 

We also introduce the useful quantity V % = A 1 + B l 
and remind that all the latin indices on the perturbed 
fields are raised and lowered thanks to the Kronecker 
flat metric 6ij. 

In the following we use the unperturbed results 

K = Mr 2 H 2 a, 
J° = 6c 2 H , J i j = 2Ua5), 
I^ = 6c 2 H6% , 1^ = 0, 7§ = -2-HaS^Sij. (B3) 

The Einstein Equation ([5]) without matter introduces the 
stress tensor 

f a p = - V CT + f a/3 
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where 



We decompose these equations into Fourier components: 



-Cs - f k 



As intermediate results wc have 

5K af3 lS = Cl (6g^g yS + g aP Sg lS ) 



8Y 0i = SYio = ci (HVt + Vl) , 

1 

2 
1 



5T l} = -M 2 Fa 2 h rjl 



sA = 



2Uc s Vi 2 Cl V! 



SJ Q l = -5J° i + —(c 1 + c 3 )B i , 

SJ,° = SJ j 8ij + — (ci + c s )B h 
a a a 



and 



<yjv = 2 Cl 
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(B5) 
(B6) 
(B7) 



ST 0l = -a BiM F — Xa Vi — FkSY^, (B8) 



(B9) 

(BIO) 
(Bll) 



(B12) 



We then obtain, 

SISi = F K [2Uc 3 V i -2c 1 Vl -2HaB i \5 Q T (B13) 
+2F K ( Cl - c 3 )d [k V t] 

with a = c\ + 3c2 + c 3 , and 

Sli Q = F#[-2ciV*' - 2c 1 HV i + 2Hc 3 V i 

+2c 3 V i ' + 6c 2 HA 1 }. (B14) 



We also get 



51% = -5%F K 5J {ij) . 



(B15) 



To obtain the perturbed Einstein equation, we then plug 
these results into the relations 

5V a I^ = do «) + d k (M*) + 2U5I& 

-{B> + UB>)I% - m 2 C2BiF K - r^o: 

sv a i° = do (6i° ) + 2U8i% - - us ki si^ 

-ST k 0j I° k - H8 ki 6I&. 



X%H>) = J2H X(t,t)Yf m \ 
T«(i,l#) = J2 J2 T(t,t)Y^ i±m) . (B16) 

to=0,1,2 ^ 

where the orthonormal modes 

y(0) j y(±l) j y(±2) (B17) 

are eigenmodes of the Laplace-Beltrami operator: 
Ay/" 1 ^ = — fc 2 Yj , I being an arbitrary set of Lorentz 
indices. For more information on these functions, see 
[13] ■ The perturbed vector field equation can then be 
written 

2\aV = F KK K' (2Hc 3 V -2c 1 V')-2c 1 F K V" 

-4F K ciHV' + V \2^-c 3 F K + U 2 {2c 3 + 2 Cl )F K 

-F K (2 Cl k 2 A + (ci - c 3 )k 2 B) + 3F K c 2 k 2 E'. (B18) 
The spin— 2 part of the ij Einstein equations gives 
k 2 E 



E + 2UE 



2F KK K'( Cl +c 3 )E = 0, 



l + F K ( Cl +c 3 ) 

(B19) 

while the spin— 1 part of the ? Einstein equations gives 
the same equation as (|B18|) and the spin— 1 component 
of the ij Einstein equation is, in space conventions, 
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2(2W + H 2 ) + M 2 Ta 2 + 4a(HF K )' + 8H 2 aF 
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+ 



1 + 2F K c 13 

ci 3 F KK K' + 2c l3 HF K 
1 + 2F K c 13 



d (i A 3l- 



d (l E 3) 



9 (i E 3) 



(B20) 



In order to compare with the results of [17| let us consider 
the gauge Ej = 0. The spin— 1 part of the ij Einstein 
equations (|E>20|) implies 



Bi — ci 3 FxAi. 
After that we can rewrite (|B18|) as 



(B21) 
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1 + F K c 13 
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and obtaining the exact expression of A in the case 
F{K) = K, 



X = 3c 
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implies in this latter case 
Cl {2UA' + A") = -- 



ci - c 3 
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A. (B23) 



It is interesting to find that in the limit where J^{K) = 
K we recover results from [l7| for (|B19IB22|) and (|B21I) . 

A quick look at (|B19I) and (|B22[) implies the existence 
of two different speeds 
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1 + F K c 13 
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(B24) 



being the sound speeds respectively of the tensor pertur- 
bation Eij and the vector perturbation Aj. These speeds 
must be positive, so that it implies constraints on the 
parameters of the theory and the function F(K). In par- 
ticular, Fk 7^ const, will imply non-trivial constraints 
between the fields of the theory and the parameters c,: 
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